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Abstract 

A new analytic result in acoustics called ‘‘Formulation 
IB, 5 ’ proposed by Farassat, is used to compute the load- 
ing noise from an unsteady surface pressure distribution 
on a thin airfoil in the time domain. This formulation is 
a new solution of the Ffowes Williams-Hawkings equation 
with the loading source term. The formulation contains a 
far-ficld surface integral that depends on the time deriva- 
tive and the surface gradient of the pressure on the airfoil, 
as well as a contour integral on the boundary of the air- 
foil surface. As a first test ease, the new formulation is 
used to compute the noise radiated from a flat plate, mov- 
ing through a sinusoidal gust of constant frequency. The 
unsteady surface pressure for this test ease is specified 
analytically from a result that is based on linear airfoil the- 
ory. This test ease is used to examine the velocity scaling 
properties of Formulation IB, and to demonstrate its equiv- 
alence to Formulation 1A, of Farassat. The new acoustic 
formulation, again with an analytic surface pressure, is 
then used to predict broadband noise radiated from an air- 
foil immersed in homogeneous turbulence. The results arc 
compared with experimental data previously reported by 
Paterson and Amict. Good agreement between predictions 
and measurements is obtained. The predicted results also 
agree very well with those of Paterson and Amict, who 
used a frequency-domain approach. Finally, an alternative 
form of Formulation IB is described for statistical analysis 
of broadband noise. 

Nomenclature 

b = airfoil semi-span (m) 

Co = ambient sound speed (m/scc) 

/ = frequency (Hz) 

/ == geometry function for airfoil surface 

F = f observed in retarded time (Eq. 2(b)) 

E* = combination of Fresnel integrals (Eq. 9(g)) 
g = vclocity-to-prcssurc transfer function 
k = u>/U, convective wave number (m -1 ) 
k = kL c / 2, reduced frequency 
L c = airfoil chord (m) 

Ci = stream- wise integral length scale 
Li = correlation length in the “»” direction 
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M = V/co, Mach number vector 
M r = M ■ r/r Mach number in radiation direction 
M u = M ■ i> Mach number in direction of z> 
p = unsteady pressure on upper airfoil surface (Pa) 
p = sound pressure radiated to observer (Pa) 
dp/ds = directional surface pressure gradient (Pa/m) 
r = x — y , sound radiation vector (m) 

R uu = correlation function of u 
Suu = spectral density of u 
t = observer time (see) 

U = uniform frcc-strcam speed (m/scc) 
u = unsteady stream- wise velocity (m/scc) 

V = airfoil velocity vector (m/scc) 

V = local velocity vector magnitude (m/scc) 
w = unsteady upwash velocity (m/scc) 

x = [xi , X 2 , £ 3 ] T , observer position (Fig. 1) 

V = [yi,y2,0] r , surface source position (Fig. 1) 

(3 = y/l-M 2 

A = co//, acoustic wave-length (m) 

p = Mk/f3 2 

xp = directivity angle (Fig. 5) 

0 = unit, inward geodesic normal (Fig. 1) 

9 = radiation direction (Fig. 1) 

po = ambient density (kg/m 3 ) 
r = t — r/co, source time (see) 

$ = random phase variable 

tjj = 2 7r/, circular frequency (Hz) 

1. Introduction 

The development of analytical methods to predict noise 
radiated from an airflow over a rigid body has been a 
subject of extensive research within the acroacoustic com- 
munity for decades. Research in this area has, in large 
part, been motivated by the desire to incorporate the re- 
sults of acroacoustic analysis into an aerodynamic design 
methodology. The present work is similarly motivated, and 
the resulting formulation should lend itself well to an en- 
gineering design- tool suite. 

The current work is specifically focused on the calcula- 
tion of far-ficld noise that results from fluctuating pressure 
on a solid surface. The acoustic analogy 1 provides a frame- 
work for the development of methods to predict noise from 
many types of sources, including noise due to unsteady sur- 
face loading. Such noise is mathematically described by 
the loading source term, or ‘‘dipole term/’ of the Ffowes 
Williams-Hawkings (FW-H) equation. 2 Because the noise 
due to an airflow over a rigid surface is often dominated 
by dipole radiation, the acoustic formulations of interest 
in this work arc determined by solutions of the FW-H 
equation with the loading source term, i.c. neglecting the 
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thickness and quadrupolc terms. 

The solution of the FW-H equation can be written in 
many forms. A new solution, called “Formulation IB” is 
presented herein. The types of noise that can be calcu- 
lated with the proposed Formulation IB arc most types of 
broadband noise, including noise that is dominated by an 
airfoil’s leading and trailing edges. Formulation IB is the 
simplest loading-noise prediction formula known to date. 
This simplicity makes the new formulation highly suitable 
for statistical analysis of broadband noise for rotating sur- 
faces. 

In Section 2, Formulation IB is derived for the specific 
ease of a flat plate in rectilinear motion. For low Mach 
numbers and distant observers, the dominant term in this 
formulation is a far-ficld surface integral that depends on 
the time derivative and the surface gradient of the airfoil 
surface pressure. The formulation also contains a contour 
integral on the boundary of the airfoil surface that includes 
the leading and trailing edges. This line integral vanishes 
along the trailing edge if the Kutta condition is imposed. 

In Section 3, Formulation IB is used to calculate the 
noise radiated from a flat plate moving through a sinusoidal 
gust of constant frequency. The unstcad}' surface pressure 
used in this test ease is an analytical result from linearized 
airfoil theory that is taken from the work of Amict. 3 ’ 4 A 
mesh refinement study is performed to demonstrate the 
equivalence of Formulation IB with Formulation 1A,° a 
previously developed acoustic formulation that is also a 
solution of the FW-H equation. Results from this analyt- 
ical test ease arc also used to examine the velocity scaling 
properties of Formulation IB, which arc found to be con- 
sistent with the results of Curie 0 and Ffowes Williams and 
Hall. 7 The directivity of the noise induced by a periodic 
gust is also examined. 

In Section 4, the single-frequency surface pressure in Sec- 
tion 3 is extended by spectral representation to serve as an 
analytic broadband source model for incident turbulence 
noise. This surface pressure is used as input to Formula- 
tion 1B to predict broadband noise to the far field. The 
resulting calculations arc compared to experimental data 
previously reported by Paterson and Amict. 8 In Section 5, 
an alternative acoustic formulation is described for statis- 
tical analysis of broadband noise. 

2. Acoustic Formulation 

Consider a flat, finite surface moving in the plane x$ =0 
along a velocity vector V. Let f(xi,X2,t) denote a geo- 
metric function that is so defined that / = 0 on the surface 
edge, and / > 0 on the interior of the surface. Let P = V/ 
denote the unit geodesic normal which lies in the plane of 
the surface, is normal to the edge, and is directed inward 
(See Fig. 1). The velocity vector V and the plate’s geom- 
etry arc related to the coordinate axes as pictured in Fig. 
1. Note that V need not be constant in space or time. The 
only stipulation on the velocity is that the motion of the 
surface is in the same plane as the surface. 

Denote by x = [aq,a72,#3] T the position of an observer, 
and by y = [yi , y2, 0] T the position of a source point on the 
plate’s surface (Fig. 1). The unsteady perturbation pres- 
sure p(y , t ) on the surface gives rise to sound that radiates 
along r = x — y to the observer. This sound is described by 



Fig. 1 Schematic for the derivation of Formulation IB. 

p'(x, t ), the perturbation pressure that arrives at the point 
( ur t , X 2 , xs) at time t. Both x and y frames of reference arc 
considered fixed relative to the undisturbed medium. For 
some of the subtle mathematical details in the following 
derivation, see Ref. 9. 

From the FW-F1 equation, the loading noise is given by 
a solution of 

= ( la ) 

where c 0 is the ambient sound speed and h is the unit 
surface normal which, for the present ease, is equivalent to 
c 3 , the unit vector in the direction of the ^3-axis. H is 
the Heaviside step function and S the Dirac delta function. 
Evaluating the divergence in Eq. 1(a) yields 

^2 r ~ V V = -p(xi,x 2 ,t)H(f)S'(x 3 ) (lb) 

where S'(xs) denotes differentiation with respect to x 3 . 

Eq. 1(b) is the wave equation with a source term, and 
its formal solution in an unbounded domain requires the 
Green’s function S(g)/4irr, where g = r — t + rjco , and t 
and r arc the observer and source times, respectively. The 
solution of Eq. 1(b) can then be written in the form 

4 Tf p(jc, t) = — f f ^T(yi, t/2, r) H(f) S'(y 3 ) dy dr 
J — oo J 5ft 3 ^ 

Now, let r — y g and integrate with respect to g. The result 
can be written 

4np(x, t) = - f -[p] nt H(F)S'(y 3 )dy (2a) 

where the subscript “ret” denotes evaluation at retarded 
time r = t — r/co, and F is 

F(y 1 : */2 ; x,t) = /(yi, y 2 , t - r/c 0 ) = [/] re t (2b) 

Integration with respect to y 3 on the right-hand side of Eq. 
2(a) yields 

4 7 Tp'(z,t)= - j ( [ p ] ret H(F) | dy 1 dy 2 (3a) 

J 5ft 2 y 3 l r J y 3 — o 
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Performing the differentiation in the integrand of Eq. 3(a) 
yields 


where p is the time derivative of pressure evaluated relative 
to an observer that is fixed with respect to the medium at 
rest, and r 3 is the third component of the unit radiation 
vector v — r/r. Clearly, then, £3 =63 ■ r = cos 9 , where 
0 is as shown in Fig. 1. The first and second terms on the 
right-hand side of Eq. 3(b) arc of the similar form Q H(F). 
In Ref. 10, it is shown that the integration of these two 
terms in Eq. 3(a) can be written 

f QH(F)d yi dy 2 = [ QdS = f _ Q dS 

Jst 2 J F > 0 «/ />0 L 1 M r Jret 

(4a) 


= [?3 p] r etH(F) 

co r 


+ — [? 3 p],etH(F) 


(3b) 


Co r 


df „ 

s7 r3p 


S(F) 


where dS is the element of the surface area of the acoustic 
planform of / > 0 . Also, M r = M ■ r is the Mach number 
in the radiation direction, where A? = C/co is the local 
Mach number vector of the surface. 

The integrated value of the third term on the right-hand 
side of Eq. 3(b) is determined as follows. This integral is 
of the form 

I = / q{y\,y2)S{F)dy l dy 2 (4b) 

The differential surface element dy\dy 2 can be written 10 


dy i dy 2 = dC dAf = ff) (4c) 

| V 2 F| 

where dC and dAf arc differential elements of arclcngth that 
arc, respectively, parallel and normal to the the surface 
edge defined by F = 0, as shown in Fig. 2. The notation V 2 
denotes the surface gradient in the y i y 2 -planc. Moreover, 
it can be shown . 10 that 


dC _ dl 
|V 2 F| ~~ [1- A/r]ret 


(4d) 


where df is an element of arclcngth along the surface edge 
defined by / = 0 . Eq. 4(b) can now be written 



Fig. 2 Differential surface element in Eq. 4(c). 


where the subscript “/ = 0 ” denotes the reference frame of 
the moving surface. Furthermore, when referenced to the 
moving surface, Df/Dr must be zero, and it follows from 
Eq. 5(a) that df/dr = — V • £>. 

All three terms on the right-hand side of Eq. 3(b) 
arc now integrated in Eq. 3(a) over the physical surface 
f(xi,x 2 ,t) > 0 , using 4(a) and 4(c). Before writing the 
final solution to Eq. 1(a), note that p is referenced to the 
medium at rest, c.g. as measured by a transducer that re- 
mains stationary as the surface passes by it. The quantity 
p can be related to dp/ dr, the time derivative of pressure in 
the reference frame of the moving surface, c.g. as measured 
by a transducer attached to the surface. This relation is 


dp - dp . . dp 

p=-f + V'-Vp=-f-V/ 
dr dr ds 


(5b) 


where dp/ds is the gradient of p in the direction of V , and 
V' is the local magnitude of V. The minus sign in Eq. 5(b) 
results from the fact that the surface gradient and velocity 
arc measured from opposite directions (Fig. 1 ). 

Incorporating all of the above results into Eqs. 3(a,b), 
the solution of Eq. 1(a) can now be written. The result is 
Formulation IB, 


47rp'(x, t) = 


+ 



( dp/ dr — V dp/ds ) cos# 
Co r ( 1 - Mr ) 


p cos 9 
r 2 ( 1 - Mr ) 


dS 

ret 


M„p cos# 
r (l- M r ) 


dl, 

ret 


(6) 


1 



q(yi,y 2 ) S(F) 


dCdF 

|V 2 F| 




to*! dC 


|Y 2 F| 

g(yi,#2) 

[l -Mr ]ret 


d( 


where = M ■ £>, the Mach number in the direction of v. 

The first and third integrals in Eq. 6 represent the sound 
radiated to the far field, whereas the second integral rep- 
resents radiation to the near field. It is noteworthy to 
consider the relative contributions of the terms in Eq. 6 , 
under the conditions of low Mach number and an observer 
in the acoustic far- field, i.c. 


Note that the surface time derivative df/dr in Eq. 3(b), 
and contained in q(yi,y 2 ) in Eq. 4(c), is referenced to the 
undisturbed medium. However, df/dr can be related to 
the material derivative Df /Dr in the reference frame of 
the moving surface, by 


Dl 

Dt 


f = o 


= g + E.V/=g + W£ (5a) 


M <1 , r > A (7) 

where A is a typical acoustic wavelength of interest. With 
respect, to M and r, the surface far-ficld integral, i.c. the 
first integral in Eq. 6 , is proportional to 1/r, whereas the 
second and third integrals arc proportional to 1 /r 2 and 
Af/r, respectively. Therefore, the far-ficld surface integral 
dominates the signal under the conditions in Eq. 7. 
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Note, again, that Eq. 6 is valid for the ease of non- 
uniform flow. Therefore, Eq. 6 can be used, as is, to 
predict loading noise from rotating surfaces. Its predeces- 
sor, Formulation 1A,° is significantly more complicated in 
its rotational form, and cannot be approximated by only 
one surface integral in the far field. Such a significant sim- 
plification for far-ficld calculations makes Formulation IB 
more suitable for statistical analysis of broadband sources 
for rotating surfaces. A statistical formulation based on 
Eq. 6 will be addressed in Section 5. However, the focus 
of the current work is the time-domain application of For- 
mulation IB , as will be demonstrated in the following two 
sections. 

3. Sinusoidal Gust of Constant Frequency 

Any noise prediction made with Eq. 6 will be only as 
good as the input surface pressure p(y. t). The current 
thinking is that such time-dependent pressure data would 
result from experimental measurement or a computational 
fluid dynamics (CFD) calculation. However, in this initial 
work, an analytic expression is used for p(y, t) to serve as 
input data for the new acoustic formulation. To this end, 
an analytic formulation from thin-airfoil theory will be used 
to describe the unsteady surface pressure that results from 
a sinusoidal gust of constant frequency. This particular 
surface pressure formulation is chosen as an example that 
will be used in this section to establish the equivalence of 
Formulation IB and Formulation 1 A° and to examine the 
velocity scaling properties of Formulation IB. In addition, 
the frequency- dependent directivity of the far-ficld sound 
produced by this sinusoidal gust is discussed. 

3.1 Surface Pressure from Thin Airfoil Theory 

Consider a rectangular flat plate, in rectilinear motion, 
as in Fig. 3. The velocity vector V = [— U, 0 , 0] T , where U 
is a constant subsonic speed. For the following examples, 
the plate’s surface and its boundary, / > 0 , arc defined by 
{0 < xi < L c ) x {— b < X 2 < b] . This surface will be 
presumed to have an unsteady pressure distribution that 
is analytically prescribed from linearized airfoil theory, as 
discussed below. 

In Refs. 3 and 4, Amict presents closed-form expressions 
for the unsteady pressure on the surface of an infinite-span, 
thin airfoil. The airfoil is presumed to move rcctilincarly 



through a sinusoidal gust. Analytical methods arc used to 
solve the two-dimensional, time- dependent linear potential 
equation by representing the solution as a product of spa- 
tial and temporal solutions. The solution is represented 
as a truncated series in which higher-order terms arc ne- 
glected (See Refs. 3 and 4 for details.). 

A complex-valued representation for the airfoil surface 
pressure is assumed to arise from a stationary gust in one 
spatial dimension. This gust can be written in the “sta- 
tionary” variable x\ — Ut as 

w{x 1 ~Ut) = wo c- ,k{ - Xl ~ ut) (8) 

where k = u/U is the stream- wise convective wave number, 
and too is the gust amplitude. This gust and the airfoil 
surface pressure that it generates arc, for now, considered 
as functions of a single amplitude and frequency. 

The unsteady surface pressure that arises due to the 
incidence of a gust of the form in Eq. 8 can be written 

A P(xi,t) = poUwo g(xi, k) c %kut (9a) 

where po is the ambient density, k = kL c / 2 is the re- 
duced frequency (based on the semi-chord), and g(xi,k) 
is a transfer function whose form is dependent on the fre- 
quency of interest. In Ref. 8, the suggested parameter to 
delineate between the low and high frequency regimes is 
p = Mk/ji 2 , where j3 = \/\ — M 2 . 

For low frequencies, p < 0.4, the transfer function is 

g(x\, k) = - l) G s (F)c'* , p < 0.4 , 

(9b) 


where k* = k/j3 2 , Gs is the classical Scars function , 11 
which, for the present work, is approximated by 


Gs(k-) 


1 

1 + 2Ak * 


+ 2nk * 


(9c) 


as suggested in Ref. 8 , and 

q(x u M) = M 2 (2 Xl /L c -l) +(1— /3)lnM+ / 31n(l+/3) — In 

(9d) 


For high frequencies, fx > 0.4, the transfer function 
is the sum of a leading-edge solution and a trailing-edge 
correction , 4 i.c. 

g(xi,k) = (91 +g2 ) c -*Mi-M) Xl/ L c+ ,/4] _ yu > 0.4 

(9c) 


where 


gi(xi,k) 


g2(xi,k) 


and 


1 

[2jrfca;i(l + M)/L c ]i 

(Of) 

-1 + ( 1 + 8 ) £*[4,t<(l - xi/L 0 )] 

['25rfc(l + M) ] i 


Fig. 3 Schematic for the constant-frequency loading noise 
problem in Section 3. 


£*«) = 


0 ( 2 7TU ) ^ 


du 


( 9 g ) 
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The quantities C { ( ) and 6’ ( ( ) arc the Fresnel cosine and 
sine integrals, and will be evaluated numerically by the 
formulas derived by Boersma. 12 The final representation 
for the unsteady surface pressure p(x\,t) assumed to be a 
real quantity, is 


p(x\, t) = 9te{ AP(a;i , t)} (Oh) 

Note that Eqs. 9(a-h) represent the pressure distribution 
on the upper surface of the airfoil and that this pressure 
is assumed to be antisymmetric between the upper and 
lower surfaces. Note, also, some differences between Eqs. 
9(c-f) and their counterparts in Refs. 3 and 8. Such dif- 
ferences include the choice of coordinate-axes origin and 
the spatial normalization employed by the author. An- 
other difference arises from the use of complex conjugates 
here, which serves to make the present notation more con- 
sistent with the eventual broadband representation in Ref. 
8, which is also employed in the Section 4. The above sur- 
face pressure was originally proposed in order to derive an 
expression for unsteady lift which was ultimately incorpo- 
rated into a frequency-domain acoustic formulation. 3,8,13 
However, in the present work, the unsteady pressure itself 
will be used as input to Eq. 6 for a time-domain prediction. 


3.2 Grid Refinement Study 

The surface pressure in Eqs. 9(a-h) is now used to numer- 
ically demonstrate the equivalence of Formulation IB and 
and the loading-noise terms of Formulation 1A J . Formula- 
tion 1A forms the basis of WOP-WOP, a rotor noise predic- 
tion code developed at NASA Langley Research Center. 14 
For the prediction of loading noise from an airfoil in uni- 
form rectilinear motion, Formulation 1A simplifies to 


47rp'{x,t) 


■ 


> f 1>0 


>f 1 >0 


dp /8 t cos 9 


cor (l - M r ) 2 J 

p cos 6 


dS 


dS 


1>0 t 


r 2 (1 - Mr ) 2 J 
P COS 0 [Mr — M 2 ) 1 


( 10 ) 


r 2 ( 1 - Mr ) 3 


dS 

J ret 


Note that the entire formulation here is integrated on the 
surface interior. At first glance, the form of Eq. 10 appears 
no more complex than Eq. 6, but only because of the sim- 
plicity of uniform, rectilinear motion. Eq. 10, as written 
above, is not applicable to a rotating surface, but Eq. 6 is. 
The loading-noise terms of the full Formulation 1A arc, in- 
deed, applicable to rotational flow, but the full formulation 
is more complex than Eq. 10. 

The far-ficld noise radiated from a thin airfoil in a one- 
dimensional, single-frequency sinusoidal gust is now calcu- 
lated, using Eq. 6 and Eq. 10. Let p \b and Pi a denote 
the sound calculated by Eqs. 6 and 10, respectively. If the 
input surface pressure p(y,t) is known analytically at any 
point on the airfoil surface, then the onty non-machinc-zcro 
error made in the numerical solution of Eqs. 6 and 10 is the 
error associated w r ith the quadrature formula that is cho- 
sen to perform the surface and contour integrations. In this 
ease, the mid-point rule is the quadrature of choice. There- 
fore, the equivalence of Eqs. 6 and 10 is demonstrated if the 
difference | p\s ~ Pi a I diminishes in mesh refinement like 


the cummulativc error expected from the mid-point rule, 
i.c. that the error is 0( A# 2 ). 

The plate’s rectangular dimensions arc determined by a 
chord length of L c = 0.5 meter and a span of 2b = 2.0 
meters. The plate is moving at a Mach number of 0.2, 
and the sound speed is taken to be 343 m/scc. The am- 
bient density po = 1.23 kg/m°, and the up wash amplitude 
is wo = 0.05 U, i.c. five percent of the free stream. The 
observer position for this test ease is x = [—1, 0, 1] , in 
meters. Fig. 3 roughly depicts this relative observer posi- 
tion, although not to scale. 

The calculation is performed for one time period of the 
surface pressure fluctuation at frequencies of / = 25 Hz 
and / = 1 kHz, with 32 time-steps in each period. These 
choices of frequency, at the prescribed observer location, 
will test both the near-field and far-ficld equivalence of the 
two formulations. Note that the transfer functions in Eqs. 
9(b,f) arc singular at x\ = 0, and the spatial derivative of 
Eq. 9(f) is singular at x\ = L c . Although both singulari- 
ties arc intcgrablc, they would cause the quadrature error 
to deviate from that of the mid-point rule which, by its 
definition, requires sufficient smoothness throughout the 
interval of integration. Therefore, the domain of stream- 
wise integration is restricted to an interval of the form 

eL c < x\ < (1 — e) L c ( 11 ) 

where c is a small, positive parameter. 

Each calculation is performed on a sequence of six sur- 
face grids: {10x40}, {20x80}, {40x160}, {80x320}, {160x 
640}, and {320 x 1280}. Grid clustering is performed near 
the leading and trailing edges of the plate in order to acco- 
madatc the parameter £ = 0.02 on the coarser meshes. The 
maximum values of | p\ B — p lA | during each time period, 
arc shown as a function of the number of grid points on a 
log-log plot in Fig. 4. The abscissa N x is the number of 
surface elements in x\. The dashed line represents a ficti- 
tious quantity whose values arc specifically calculated to be 
directly proportional to N~ 2 . Clearly, the slopes of both 
calculations arc visibly parallel to a slope of 2, thereby 



Fig. 4 Grid refinement validation for equivalence of For- 
mulations 1A and IB. 
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demonstrating that 

\p'ib - p'ia I = 0(N~ 2 ) = O(Axi ) , 

which is the cummulativc error expected from the mid- 
point rule of integration. 

3.3 Velocit y Scalin g Laws 

Attention is now turned to the way in which the intensity 
of the far-ficld noise, as predicted by Eqn. 6, will scale as a 
function of velocity, when the surface pressure is described 
by Eqs. 9(a-h). The physical dimensions of the plate for 
this exercise arc the same as for the test ease above. The 
conditions of a far-ficld observer in a low-Mach-numbcr 
flow (Eq. 7) will be assumed. The Mach-numbcr range of 
interest is 0.01 < M < 0.2. In addition, the proportionality 
of the acoustic intensity to velocity will be determined ac- 
cording to two different conditions placed on the frequency 
/. These two conditions will delineate between compact 
and non-compact sources. 

First, it will be assumed that the source is compact, i.c. 
L c A. This condition will be achieved by requiring the 
frequency to be proportional to velocity, / ~ U/L c , for a 
sufficiently low range of frequency. For 0.01 < M < 0.2, 
/ = 80 Hz is referenced to M = 0.2, and / is then made 
a function of velocity while requiring a constant Strouhal 
number f L c /U « 0.583. The frequency range is, then, 
AH z < / < 80 Hz. Note that p < 0.4 throughout this 
range. 

The plate’s physical dimensions arc as in the above mesh 
refinement problem. The calculation is performed on a 
100 x400 surface grid, with the stream- wise integration in- 
terval restricted as in Eq. 11, with sufficient grid clustering 
near the leading and trailing edges to allow for e = 0.003. 
The observer is chosen at a distance of 100 meters, directly 
above the plate’s center, i.c. x = [0.25,0, 100 ] 1 in meters. 
This location places the observer in the acoustic far-ficld 
for the entire range of frequency. 

A separate calculation is run for each Mach number 
and its corresponding frequency. The upwash amplitude is 
we = 0.05 U for each of 50 equally spaced Mach numbers 
between 0.01 and 0.2. The surface pressure in Eqs. 9(a- 
d,h) is used as input to equation to Eq. 6 to predict the 
far-ficld sound p'(x,t). Each calculation is performed for 
one acoustic period T of the corresponding frequency, with 
64 time-steps. The average intensity I{x) of the acoustic 
signal at the observer x, assuming spherical spreading, is 
then calculated by 

I f T Lp'OM )] 2 dt 

T Jo Po co 

The average acoustic intensities for a compact source, as 
a function of Mach number, arc represented as squares in 
Fig. 5. The slope of these results on a log-log plot can be 
visually determined by proximity to the dotted line whose 
slope is precisely six at every point. This U G proportional- 
ity is consistent with Curie’s 4 result, as expected from the 
conditions placed upon the calculations. 

The demonstration of a velocity scaling law is now de- 
sired for a non-compact source i.c. for /U > 0.4. Therefore, 
the restriction that / ~ U / L c must be lifted, so that / is 


I(x) 



Fig. 5 Velocity scaling properties of Formulation IB. 


independent of U . The simplest such condition is that / 
is constant, in which ease the Strouhal number remains a 
function of U . A series of calculations is again performed, 
as above, with the only parameter change being that the 
frequency is held constant at 1 kFlz throughout the range 
of Mach number. The input surface pressure for this ease is 
given by Eqs. 9(a,c-h). The computed acoustic intensities 
at 100 meters arc represented by circles in Fig. 5. In this 
ease, the acoustic intensity scales approximately as U° , a 
result that is consistent with Ffowes Williams and HalF . 

3.4 Directivity 

As a final exercise in this section, the directivity of a 
single-frequency source is examined. The radiated noise 
p'(x , f ) is calculated at many locations on a circular arc in 
the plane X 2 = 0 that is centered on the geometric center of 
the plate’s upper surface, as shown in Fig. 6. The arc tra- 
jectory ( r ■ «/’) is determined by r = 3 meters and 0 < ip < tt. 
The directivity is determined by the peak pressure ampli- 
tude Up'll that, is calculated at each position on the circular 
arc, during one period in time for a given frequency. The 
flat plate’s dimensions and surface discretization arc as in 



Fig. 6 Schematic for directivity calculation. Observer on 
circular path in plane a; 2 = 0. 
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\p' Imax, Pa. 

Fig. 7 Directivity for two sources of constant frequency. 
2-kHz results arc sealed by a factor of 2.5. 


upstream direction. Far-ficld noise measurements of the 
incident turbulence on the airfoil arc determined by sub- 
tracting microphone measurements, with and without the 
model, at each of five tunnel speeds: 40, 60, 90, 120, and 
165 m/s. 


4.2 Broadband Analysis 


For prediction purposes, the airfoil is modeled as a flat 
plate in a periodic gust that gives rise to an unsteady 
surface pressure that is a broadband extension of the an- 
alytic formulation in Section 3. The airfoil geometry is 
oriented with respect, to the coordinate axes as in Fig. 3, 
with { 0 < aq < L c } x {— 6 < x 2 < 6), where L c = 0.23 
m and 26 = 0.53 m. As encountered by the airfoil surface 
in the x'i-x '2 plane, the normal component, of the turbulent 
velocity field can be written 


the previous example, and there arc 128 time-steps in a pe- 
riod. The observer path, 0 < */> < n, is discretized into 128 
equally spaced locations. The frcc-strcam Mach number is 
0.2 and the gust amplitude is wo = 0.05 U. Fig. 7 shows 
the results, in polar form, for frequencies of 1 kHz and 2 
kHz. The higher-frequency results axe sealed in order to 
visualize both loci on the same plot. As expected, there arc 
twice as many lobes in the 2-kHz solution. Also, note the 
frequency-dependent positions of the lobes with respect to 
a fixed observer. For example, an observer at </> = jt/ 2 
will receive a signal that is near the peak amplitude for the 
central lobe of the 1-kHz signal, whereas the 2-kHz sig- 
nal is near a local minimum for the same observer. This 
frequency-dependent characteristic of directivity is men- 
tioned here for future reference in the following section. 


4. Broadband Prediction with 
Com p arison to Experiment 

The analytic surface pressure in the previous section is 
extended to model a broadband source on a slender airfoil 


at zero angle of attack. This broadband surface pressure 
is used as input to Formulation IB to predict far-ficld 
radiation, and the results arc compared t.o experimental 
measurements. Note that. Amict 15 has previously proposed 
a broadband solution to this problem in the time domain, 
using the transfer functions in Eqs. 9(c-g). His resulting 
solution was a Fourier transform of the frequency-domain 
solution in Refs. 8 and 13. In the current, work, the high- 
frequency formulation in Section 3 is explicitly extended to 
a broadband source application with user-specified spectral 
content. 


v(x\,X 2 ,t) = 


/ oo pc 

-oo j — l 


'j(ki,k 2 ) i 


-c/t) +^] dkl dk 2 


where w[ki,k 2 ) is the gust, amplitude wave-number 
spectrum, defined by the inverse Fourier transform of 
w(xi,x 2 ,t). The complex-valued, unsteady surface pres- 
sure arising from the incidence of a turbulent, velocity field 
of this form is given by 

AP(:tl,I 2 ,f)= (12) 

poo p CO 

poU / w(k\,k 2 )g(x\ ,ki ,k 2 )c' kl 1 ‘~ k2X2 ^ dk\dk 2 


The expression for surface pressure in Eq. 13 is simpli- 
fied by the following reasoning. Amict argues in Ref. 13 
that, for an observer in the plane x 2 = 0, the only spanwisc 
wavenumber that contributes significantly to the far-ficld 
sound is k 2 = 0. His conclusion is derived mathematically, 
in the frequency domain, for the limiting ease of an airfoil of 
infinite span. From a physical standpoint, this conclusion 
makes sense for an observer in a location that is symmetric 
to the airfoil span. The effect on the far-ficld acoustics of 
any gust that is skewed to the airfoil leading edge by some 
angle a will be canceled by another gust skewed at an an- 
gle of —a. In the ease of a finite-span arifoil, Amict argues 
that this simplification is still valid as the quantity Mkib 
becomes large, i.c. in the high-frequency limit.. Follow- 
ing this line of reasoning, the x 2 dependence in A P also 
vanishes, by Eq. 12, and the surface pressure becomes 

/ OO 

w(ki ,0) g(xi ,0) c l ^ kl L 1 ^ dk\ (13) 


4.1 Ex p eriment - Incident Turbulence Noise 

The experiment that is modeled in this section is re- 
ported by Paterson and Amict in Ref. 8. A NACA 0012 
airfoil is placed between two vertical plates, at zero angle 
of attack, in the test section of an open-jet wind tunnel. 
The airfoil has a chord length of 0.23 m and a span of 
0.53 m. Turbulence is generated by a grid upstream of the 
airfoil. Noise propagates from the test section into an anc- 
choic chamber that is instrumented with six microphones. 
The microphones arc located on the tunnel centerline, on 
an arc of raidus 2.25 m, relative to the airfoil’s geometric 
center. The microphone locations on this arc arc at angles 
of 70, 90, 105, 120, 130, and 140 degrees, relative to the 


The evaluation of the surface pressure in Eq. 13 is accom- 
plished by first recognizing the turbulent fluctuations as a 
stochastic process. This process can be approximated by a 
truncated series whose limit exhibits the required relation- 
ship between the autocorrelation and the power spectrum 
of that process. (See, for example, Ref. 16. ). This rela- 
tionship is achieved by evaluating the spectral amplitudes 
w(k\, /t 2 ) as a function of the power spectral density (PSD) 
of w(x\ , X 2 , t). To this end, the infinite wave-number do- 
main, — oo</ci<oo, in Eq. 13 is integrally discretized and 
truncated such that k i,— iv < k \, n < • The largest wave 

number k\ y N represents an ; ‘uppcr-cut-off” wave number, 
beyond which the spectral density amplitude is considered 
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negligible or is out of range of experimental measurement. 
The unsteady surface pressure in Eq. 13 is then approxi- 
mated by 


N 

APfii, t) w ppU y, A n ,a c 1 ' 5 ’’* g(xi, fci,n, 0)c ! ’ 

-N 

ki, n = n Aki , n = 0, =tl, ±2, . . . , ±iV 
A/ci = /iV 

The phase angles {4> 7 z} arc independent random variables 
uniformly distributed on [0, 2tt\. The transfer function 
<7(37 1, ki }U , 0) is an extension of the high-frequency function 
in Eqs. 9(c-g), and takes the slightly modified form 

g{x l,fcl,„,0) = (31 +g2 ) c -‘[ 2 ^U-M)«/ic + ^/4-S 1 ,„] 


where gi and <72, now functions of /ci >nj arc otherwise iden- 
tical to Eq. 9(f), 




92 ( 3 ;i, fcl,„) 


1 

[2xAi,„aji(l + M)/L c ]i 

(14b) 

-1 + ( 1 + i ) £*[4ju n (l - xi/Lc)] 
[2n-fci,„(l-|-Af)]2 


and £* is the same complex combination of Fresnel inte- 
grals as in Eq. 9(g). The low-frequency transfer function 
is not used in these broadband predictions because the ex- 
perimental facility is ancchoic for frequencies above 200 Hz, 
and the parameter p is greater than 0.4 at this frequency 
or above, for all five tunnel speeds. 

The spectral coefficients {A„o ] arc evaluated by 


A n ,0 — [ Sww (k\ } n , 0) A &1 A&2 ] 2 ( 14c ) 


where S ww (k u k 2 ) is the two-component power spectral 
density of to. As this grid-generated turbulence is assumed 
to be homogeneous and isotropic, the PSD is evaluated by 
the von Karman formula. 1 ' The formula for S mw (ki, kz) 
that is used in the present calculations is derived from von 
Karman’s energy spectrum in Appendix I of Ref. 13. 


S ww (k i , /C2) 


4 u 2 k\ -f fcf 

Q* k l [1 +kj + kl]i 


(14d) 


where u 2 is the stream- wise turbulent energy, and 


ki = 


kj_ 

k e 


k = v^£l|) 

' A r (l) ’ 


and C\ is the stream- wise integral length scale 


(14c) 


£1 


R WW (# 1 ) dx 1 


then set to zero for the present calculations. Values for 
u 2 and Ci arc determined by measurement. In Ref. 18, 
Fink reports that the turbulence intensity that results from 
the grid in question can be approximated by the empirical 
formula 



0.04 


U 

Ure f 


(14f) 


where the reference speed is U re f = 60 m/s. Fink 18 also 
reports a measured value for the integral length scale as 
Ci = 3.175 cm 

All of the above expressions and measurements arc in- 
corporated into Eq. 14(a). The final representation for the 
unsteady broadband pressure on the airfoil’s upper surface 
is then given by the real part of Eq. 14(a). This broad- 
band surface pressure is used as input to Formulation IB 
to predict the far- field noise p'(x, t). 

Using symmetry arguments and algebraic manipulation, 
the indicial bounds for the surface pressure’s spectral rep- 
resentation arc altered so that the domain includes only 
positive wavenumbers. The resulting real-valued surface 
pressure on the airfoil’s upper surface can be written 


p(xi, t) = 2p 0 U ^2 An,o[Bn cos(fci,„t/f + 5>„ - . 


1 

+ D n sin(fci,„(7 1 + $ n - a n ) ] 


(15a) 


where the upwash amplitudes A„,o arc evaluated by Eq. 
14(c) and S ww (ki,k 2 ) is described in Eqs. 14(d-f). The 
quantities B„ . D n , and a n arc given by 


B n 

D n 

& n 


g\(00l,kl t n) + 


c(6Q + s(t«) - i 


[2wfci,„(l + M)] 2 
S(in) - C(( n ) 


[2jrfei,„(l + M)] 2 


2 p n ( 1 + M)— k\ >n -(- — 

L/c 4 


(15b) 


where <7i(a?i,&i,n) is evaluated in Eq. 14(b), and C(£„) and 
«S(^ n ) arc the Fresnel cosine and sine integrals in Eq. 9(g) 
with = 4 fj n ( 1 — xi/Lc )• The summation in Eq. 15(a) 
begins at n = 1 because S ww ( 0, 0) = 0, by Eq. 14(d). 

At this point, the value of A&2 in Eq. 14(c) is unknown 
because there is no explicit span-wise integration, thereby 
giving rise to an adjustiblc constant. This constant is one 
of scale only, and found to have no effect the shape of the 
far- field spectrum. Furthermore, this scale factor is found 
to be constant for all data points, i.c. is independent of 
all parameters considered (tunnel speed, frequency range, 
band- width, airfoil span, etc.). 


4.3 Time-Domain Predictions 


where Rw w is the upwash correlation function defined bj' 

Rwio(r) = w(xi. t) w(xi + r, t) 

1 f T 

= lim — / «;(#!, t) w(x\ T r, t) dt 

T—*oo 1 J 0 

Note that the two-component formula in Eq. 14(d) was 
derived in Ref. 13 by integrating the von Karman energy 
spectrum over all &3 components. The A?2 component is 


The lower frequency bound, and therefore the funda- 
mental frequency, for all five calculations is chosen at 10 
Hz. The upper frequency for the predictions is chosen ac- 
cording to the upper frequency for which measurements arc 
available for each tunnel speed. For U = 40, 60, and 90 
m/s, the upper bound is fn = 2.5 kHz. The upper bounds 
for U = 120 and 165 m/s arc /jv = 3.5 kHz and 4.5 kHz, 
respectively. For all five calculations, the numerical band- 
width is A/ = 10 Hz. Each calculation is performed for one 
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period of the lowest frequency, i.c. T = 0.1 second. The 
numerical solution is sampled at the Nyquist frequency, i.c. 
At = T/2N. The calculation is performed on a 100 x 230 
surface grid, with the stream-wise integration interval re- 
stricted as in Eq. 11, with sufficient grid clustering near 
the leading and trailing edges to allow for e = 0.003. 

The experimental microphone position for which com- 
parisons arc made is at a distance of 2.25 m from the model, 
and at an angle of 90 degrees relative to the model’s geo- 
metric center. The coordinate system for the calculation 
is such that the a? 2 -axis is coincident with the center-span 
line, so that the microphone position is in the plane X 2 = 0, 
as in Fig. 6. The measured observer position for the pre- 
diction is. then, x = [0.115, 0, 2.25] t in meters. 

The position of the microphone relative to the airfoil is 
corrected for refraction due to the presence of a shear layer 
that forms downstream of the upper lip of the square nozzle 
exit and is positioned between the model and the micro- 
phone. This correction is based on geometrical acoustics 
with an assumption of a zero-thickness shear layer, and is 
reported in Ref. 19. Shear-layer corrections that arc based 
on such formulations 13 arc reasonable for the present ease 
with the microphone directly above the source. The re- 
quired correction in the microphone position is significant. 
At a measured angle of 90 degrees, the corrected angles 
ranged from approximately 84.5 degrees for U = 40 m/s to 
68.2 degrees for U = 165 m/s. 

In addition, the amplitude of the radiated noise is also 
corrected for the presence of the shear layer, although at, a 
measured angle of 90 degrees, the amplitude correction is 
not significant, especially for the lower tunnel speeds. The 
computed sound pressures p'(x, t ) were corrected by factors 
ranging from approximately 0.997 for U = 40 m/s to 0.942 
for U = 165 m/s. The microphone position is corrected for 
the shear-layer in a pre-processing step. After the far-ficld 
noise is calculated at the corrected position, the results arc 
then post-processed for amplitude correction. In this way, 
the corrected predictions can be compared to the experi- 



Fig. 8 Predicted far-ficld signal, U = 165 m/s. Micro- 
phone at 90°, 2.25 m above airfoil center. 


mental results “as measured.” 

Fig. 8 shows the far-ficld signal p(x , t) that is predicted 
by Formulation IB, for one fundamental period in time, at 
the experimental microphone location, for a tunnel speed 
of 165 m/s. Shear- layer corrections for amplitude and di- 
rectivity arc included in this plot. In order to compare with 
experimental measurements in Ref. 8, the time-domain re- 
sults from the numerical predictions, for all five tunnel 
speeds, arc Fourier analyzed and converted to the fre- 
quency domain. The resulting frequency-domain solution, 
p\fn )j is used to compute the sound pressure level (SPL) 
spectrum of the far-ficld radiation. These sound pressure 
levels arc determined by 

SPL(/„) = 10 log (16) 

where the reference pressure is p re f = 2 x 10~' J Pa. The 
SPL’s arc converted to a 1 Hz band-width by reducing the 
values in Eq. 16 by 101og(A/). This narrow-band conver- 
sion is consistent with the experimental SPL’s which were 
measured at a band-width of 55.7 Hz and reduced by 17.5 
dB. 8 

The predicted far-ficld spectral density for the five tun- 
nel speeds is shown Fig. 9 along with experimental mea- 
surements from Ref. 8. The solid symbols represent those 
measurements for which the difference between the noise 
with and without the airfoil model was considered too 
small, and arc therefore subject to greater uncertainty. The 
agreement with the measured data is very good. The no- 
ticeable “humps” in the predicted spectra arc, most likely, 
the result of the changing placement of lobes, as a function 



Fig. 9 Predicted and measured far-ficld noise spectra. 
Microphone at 90°, 2.25 m above airfoil center. Experi- 
mental data reproduced from Ref. 8. Solid symbols denote 
low signal-to-noisc ratio. 
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of frequency, relative to the observer. This phenomenon 
was described in Section 3 (Fig. 7). That the current pre- 
dictions arc nearly identical to the predictions in Ref. 8 is 
not surprising, because the current predictions rely upon 
the same unsteady surface pressure formulation as those in 
Ref. 8. 


5. Statistical Formulation 

Often, when acroacoustic experiments arc performed, 
surface-pressure correlations arc extremely useful in char- 
acterizing noise due to an airflow over a model. Under 
certain conditions. Formulation IB is readily transformed 
into an expression that aids in statistical analysis of broad- 
band noise. Recall that, in the ease of low Mach number 
and with an observer in the acoustic far field (Eq. 7), the 
signal is dominated by the first integral in Eq. 6. Also, 
recall that the spatial and temporal derivatives of pressure 
in this integral can be written as a single time derivative of 
pressure when evaluated in a reference frame that is fixed 
relative to the medium at rest, as in Eq. 5(b). In this fixed 
reference frame, the far-ficld sound at low Mach number is 
approximated by 


4 *p'{x,t) 


f 

p cos 9 

jf> 0 

_c 0 r (1 - M r )_ 


dS 

ret 


(17) 


In addition, assume that the observer is many correlation 
lengths into the far field, i.c. r 1;, where 



In this ease, the values of r, M r , and 9 arc nearly invariant 
within a correlation area At of size by ( 2 , and Eq. 17 
can be re-written as 


4 if p\x,t) ss 


cos 9k 

co rk (1 - Mr* ) 


/ p(y, 
J A. 


t) dS (18) 


where rk, 9k, and M,- k arc constant values chosen to replace 
their nearly invariant counterparts in retarded time, within 
each correlation area. 

If the autocorrelation operator is applied to Eq. 17, the 
result is 


I671- 2 R p i p i (r) = 167r 2 p'(x, t)p'(x, t + t) 


(18a) 




cos Ok 


co r k (1 — Mr k )\ J A 


Rppi’J', n, T ) drjdy 


where 


Rpp(y'>n, T ) =P(y,T)p{y + f),r+ r) (18b) 

The ability to measure such time-derivative correlations, 
relative to the medium at rest, is already in hand. 20 The 
potential usefulness of this alternate formulation is clear. 
After the autocorrelations in Eq. 18(b) arc determined ex- 
perimentally, and used as input to Eq. 18(a), then the 
Fourier transform of Eq. 18(a) produces the far-ficld noise 
spectrum. The further development and testing of this 
statistical formulation is a topic of ongoing research for its 
potential application to trailing edge noise prediction. 


Concluding Remarks 

A new formulation for the solution of the loading term 
of the Ffowes Williams- Hawkings equation has been de- 
rived. The potential usefulness of time-domain solutions 
for acoustic predictions with Formulation IB has been 
demonstrated. This new far-ficld formulation has some 
advantages over previous formulations. The formulation 
is both simple and has broad application, including the 
ease of non-uniform flow. In addition, the dominance of 
only one term in this formulation makes this solution much 
easier to pose in correlation-function form for statistical 
analysis of broadband noise. Such an alternative formula- 
tion can be used to aid acoustic experiments where surface 
pressure correlations arc measured. 
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